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Abstract: We develop a formalism that allows the computation of the quantum 
effective potential of a scalar order parameter in a class of holographic theories at 
finite temperature and charge density. The effective potential is a valuable tool for 
studying the ground state of the theory, symmetry breaking patterns and phase 
transitions. We derive general formulae for the effective potential and apply them to 
determine the phase transition temperature and density in the scaling region. 
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1. Introduction 



Many recent applications of the AdS/CFT correspondence (in QCD or condensed 
matter theory) entail the analysis of an effective holographic theory (EHT), namely 
a classical theory of gravity coupled to a set of matter fields, e.g. Einstein theory 
coupled to a set of real or complex scalar fields and a set of gauge fields. Such 
an approach is the core of bottom-up approaches to holographic problems. It was 
advocated in [1] that this approach is well suited to study the IR asymptotics of 
holographic theories and a host of new (generalized) classes of quantum criticality 
were found at finite density. It is also in accord with alternative setups for the 
holographic flow in the IR, [2] . 

The choice of interactions between the fields in the Lagrangian is limitless in 
bottom-up approaches, and the usual strategy involves the choice of a specific well- 
motivated subset of generic interactions parametrized by a number of arbitrary func- 
tions. For instance, the minimal set of a single real or complex scalar field with a 
f/(l) gauge field has been the standard paradigm in constructions of phenomeno- 
logical models of YM [3, 4] and s-wave superconductors at finite temperature and 
charge density, [5]. 

One of the central objects of interest in such systems is the free energy, which can 
be computed using standard methods by evaluating the on-shell gravitational action. 
This quantity, as a function of external sources, is related by a Legendre transform 
to the effective action of the theory. The effective action contains all the information 
needed to determine the vacuum structure and other generic correlation functions. 
An important part of the effective action is the (quantum) effective potential. It is an 
important ingredient in deciding spontaneous symmetry breaking and superfluidity 
(or superconductivity) at finite density. 

Solving the full non-linear set of gravitational equations at finite temperature 
T and charge density p (in generic systems of the type outlined above) is a highly 
non-trivial task that typically requires numerical work. The purpose of this paper 
is to give compact expressions for the quantum effective potential in holographic 
theories in a form that facilitates further numerical computation in concrete models. 
This is achieved by generalizing methods applied previously to zero temperature and 
zero charge density solutions [6, 7]. We introduce a properly defined superpotential 
function W and use it to reduce the second order differential equations to first order. 
This reduction simplifies the differential equations that determine the background 
solution. The quantum effective potential of the theory (which is determined from 
the on-shell gravitational action) can be expressed compactly in terms of W and the 
entropy function. We demonstrate how finite temperature and charge density effects 
enter into this expression. 

Perhaps the most appropriate way to think about the superpotential function is 
within the context of the Hamilton- Jacobi treatment of holographic renormalization 
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[6, 7]. In this framework, the on-shell gravitational action coincides with Hamilton's 
principal function Sh which obeys a first order differential equation, the so-called 
Hamilton- Jacobi equation. A solution of the Hamilton- Jacobi equation determines 
a set of first order flow equations that specify the background solution. In this 
context, the superpotential function W is the potential term in Sh, and for static 
and homogeneous solutions the first order equations are simply the Hamilton-Jacobi 
equation for W and the corresponding Hamilton-Jacobi flow equations for the fields. 

This elegant reformulation of the second order equations of motion is generic 
and applies also to finite- (T, p) solutions. The practically difficult part of the imple- 
mentation of this formalism in this case is to determine the appropriate solution of 
the Hamilton-Jacobi equation. Unfortunately, a generic ansatz for Sh that captures 
finite temperature and density solutions, even static and homogeneous ones, is not 
known and it is not obvious how to find it. For that reason, in the present paper we 
propose a different less ambitious approach. 

The approach that we adopt keeps the zero-(T, p) first order flow equations in 
terms of a function W intact. With this ansatz we examine the extent to which the 
second order equations reduce to a first order system and find that a full reduction 
to a first order system is possible for finite T, but zero p. In the more general case 
of arbitrary T and p some equations remain second order. The function W that we 
define in this manner is the standard superpotential that coincides with Hamilton's 
principal function Sh in the T = p = case. For general T, p, however, it does 
not coincide with Sh and the relation between the W that we define and Sh is less 
straightforward. 

Summary of main results and structure of the paper 

We apply the above strategy to a rather general class of Einstein-Maxwell-dilaton 
theories captured by the action (2.1). These theories are relevant in holographic 
descriptions of strange metals and other non-Fermi liquids and phenomenological 
models of the glue sector in holographic QCD. The reduction of the second order 
equations of motion is discussed in section 2 and the resulting equations are sum- 
marized in eqs. (2.21a)-(2.21c). A compact general expression of the free energy at 
a general renormalization (RG) scale M in terms of the function W is derived in 
section 3.2 (eq. (3.8)). This equation is one of the main results of the paper. 

The removal of M by taking M — )■ oo and the analysis of the corresponding 
divergences is performed in sections 5, 6 providing compact expressions for the UV 
quantum effective potential. We use this result in section 7.1 to determine the critical 
temperature Tc at p = that separates a normal phase without a scalar condensate 
from an ordered phase with a condensate. A general expression for the critical 
temperature at any charge density is presented in section 7.2. It is known [8] that 
this system can exhibit a quantum critical point, i.e. a critical point with = 0. 
We determine this point and related scaling exponents analytically. 
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The case of Einstein-abelian Higgs actions with a complex scalar field, which are 
relevant in discussions of s-wave holographic superconductors [5] is more complicated. 
We have included the pertinent formulae in appendix A. 

As a preliminary illustration of a qualitatively different set of theories we consider 
analogous computations in Einstein gravity coupled to a scalar field with a non-linear 
DBI-inspired action (see eq. (D.l)). Similar actions appear in discussions of the fiavor 
sector in models of holographic QCD [9]. A novelty in this case is the presence of 
multiple branches for the superpotential function W. The relevant discussion, which 
focuses exclusively on zero temperature, appears in appendix D. 

The detailed analysis of applications of this formalism is beyond the scope of the 
present paper. A brief discussion of potential applications appears in the concluding 
section 8. Additional useful technical details are relegated to the appendices B, C. 



2. Reduction of the equations of motion 

In this section we describe the reduction of the second order gravitational equa- 
tions for a generic Einstein-Maxwell-dilaton model. We consider solutions at finite 
temperature and charge density. 

2.1 First order equations 

Our first example is provided by the following Einstein-Maxwell-dilaton action in 
d + 1 spacetime dimensions 



Igh = 2M^-i Jg^ d'^xVhK . 



(2.1) 



This action describes the dynamics of a U{1) gauge field (with field strength 
F = dA) and a real scalar field coupled to Einstein gravity. The boundary term Igh 
is the standard Gibbons-Hawking term needed to make the variational problem well- 
defined. As such this action describes the grand canonical ensemble. The Lagrangian 
is parameterized by two functions {Z, V) of the scalar field 0. We assume that these 
functions have the following expansion around = 

oo oo 
Z(0) = l + ^Z„02n^ y^^^^J2^^^2n _ ^2.2) 

n=l n=0 

We are interested in solutions that asymptote to AdSd+i- In the asymptotic 
region where — )■ the first few coefficients of the expansion of the potential V are 

vo = d{d-l) , vi = — (2.3) 
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in units where the AdS radius is set to one. m denotes the mass of the scalar field 
(j). The standard AdS/CFT dictionary imphes that the bulk scalar field maps to a 
real scalar single-trace operator on the boundary with UV scaling dimension A such 
that 

m2 = A(rf-A). (2.4) 

For in the range — 1, ^) there are two sensible values of A that obey this 
equation, A-t (by convention A_ = d — A_)_ < A+). In all other cases, only A+ is 
allowed. 

Since we are interested in solutions with finite temperature and chemical poten- 
tial we set 

dip 

ds^ = e2^(") (-f(u)dt^ + dx'dx') + —— , A = At(u)dt , = 0(m) . (2.5) 

fw 

We are working in domain wall frame coordinates where the AdS solution has 

A{u) = -u , f{u) = 1 , At{u) = , 0(m) = . (2.6) 

In general asymptotically AdS solutions the UV AdS boundary lies at u — > — oo. 

In this frame the second order equations of motion reduce to the following set of 
differential equations 

2(rf- l)i + 02 = , (2.7b) 
/ + dAf - e-^^ZAl = , (2.7c) 

{d - l)Af + Ud - l)A^ - \A f~V + \ze-^^A\ = . (2.7d) 



We use the convention 

. _ d 

du 

The first equation (2.7a) can be solved trivially to obtain 



(2.f 



The integration constant p is the charge density for the U{1) gauge field. 
By introducing a function W such that 

= W^'(0), ':=^ (2.10) 

we notice that the second equation (2.7b) is solved by setting 

i = (2,11) 
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Equivalently, 



A{(j)) = Ao- 



1 



2{d-l] 



~ w{4>) 

am — , U = Un + 

W'{<P) 



(2.12) 



where Aq = A(0o)- 

Aq is an arbitrary value of the scale factor at an intermediate position in the 
radial direction. Its choice corresponds to the definition of an RG scale M = at 
which the scalar operator is defined. Therefore 0o is the effective "source" field at 
the RG scale associated to Aq. 

The third equation (2.7c) takes the form 



W'iW'f 



dWW 
2{d-l) 



f 



This can be integrated once to obtain 



W'f 



-dA 



D + p" 



d(/. 



4>o 



^{d-2)AzW' 



(2.13) 



(2.14) 



D is an integration constant that can be fixed in terms of the temperature T and 
the charge density p by requiring regularity of the geometry at the horizon u = Uh 



A-kT . 



(2.15) 



U=Uh 



Choosing a negative sign for / at the horizon we obtain as a monotonically in- 
creasing function of u} With this choice eq. (2.15) translates to 



D + p" 



^{d-2)AZW' 



-inTS 



where 



is a quantity proportional to the entropy. Eq. (2.14) becomes 

d0 



W'f 



-dA 



-AtxTS + 



<l>h 



^{d~2)AzW' 



Integrating once (with the normalization condition^ /(0o) = 1) we obtain 



/(0) = 1 + 



^ -dA 
W 



-AttTS + p2 



gid-2)AZW' 



Finally, the last equation (2.7d) becomes 



2{d 



W") f- WW'f = 2V 



(,2{d-l)AZ ■ 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



^In cases where we expand around an IR AdS solution it is appropriate to pick a positive sign 
for /. 

^This condition imposes that the metric at the RG scale is the standard flat Minkowski metric. 
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Summary of equations 

The ansatz (2.10) has allowed us to reduce some of the second order equations to a 
system of first order equations. A solution at finite temperature and charge density 
can be found by solving (typically numerically) the system of coupled equations 



T^W' = nW , 7^ := e 
dWW 



W'iW'f 



2{d 



r 



-2(d-l)A 



2{d-l) 



] / - WW'f = 2V-^ 



(2.21a) 
(2.21b) 

(2.21c) 



for the unknown functions 7^, /, W . The profile of the scalar field and gauge field 
component At can then be determined from eqs. (2.9), (2.10). 

2.2 Special cases and other comments 

In the special case of zero charge density, but arbitrary temperature, the second 
equation (2.21b) can be replaced by the first order differential equation (2.14) 

W'f = De-'^^ . (2.22) 

In this case we reduce the full set of equations to a set of first order equations. 

Moreover, at zero temperature and charge density / = 1 and the full solution 
can be determined by solving a single first order equation for the function W 



dW^ 



W'^ = 2V . (2.23) 



2{d-l] 

W is the standard superpotential function in this case, the potential (derivative- 
independent) term in Hamilton's principal function. In this language eq. (2.23) is 
the Hamilton- Jacobi equation for static and homogeneous configurations. 

In the more general case of non-zero temperature and charge density, the function 
W, as defined here, does not express the full potential contribution to Hamilton's 
principal function. Solving the Hamilton-Jacobi equation for static and homogeneous 
configurations at finite temperature and charge density remains an interesting open 
problem. 

Finally, a more general case of considerable interest for applications is the case 
of the Einstein-abelian Higgs model expressed by the action 



R - - hd<pr - Jmae - qA)' + v{<p) 



+ Igh . 
(2.24) 
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This action describes the dynamics of a U{1) gauge field (with field strength 
F = dA) and a complex scalar field $ := 0e*^ {6 G [0,27r)) coupled to Einstein 
gravity. The constant q denotes the U{1) charge of the scalar field $. 

For solutions with a constant profile for the angular field 6 we can consistently 



set dfj,6 = into the action to obtain 



+ Igh (2.25) 



which introduces a new set of interaction terms into the action (2.1). The equations of 
motion of this action and the first order reduction for solutions at finite temperature 
and charge density in the approach of the previous subsection are summarized in 
appendix A. 



3. General expressions for the holographic effective potential 
3.1 Introductory comments 

We are now in position to compute the effective potential. In order to fix the notation 
let us recall first a few standard facts. 

On the field theory side the free energy is defined, as a function of the sources, 
with a path integral of the form 

Z[J] = e-^[^l= f e-f"'^^^^'''^ (3.1) 



where C is the field theory Lagrangian. For concreteness, we focus here on a single 
scalar operator O and its source J. The vev of the operator O is obtained from a 
functional derivative of the free energy as follows 

(O). = ^ . (3.2) 

The effective action F, which is a function of the vev (O), is the Legendre trans- 
form of the free energy J-" 

T[{0)j] = ^[J] - / d'^a; J{x){0{x))j (3.3) 



in terms of which the source can be expressed as 

On the gravity side, and the standard 'quantization'^ of the dual scalar field 
0, the free energy J^[J] is computed by evaluating the on-shell gravitational action 



^The so-called standard 'quantization' corresponds to the situation where the scalar operator O 
has UV scaling dimension A+ in the notation below eq. (2.4). In the alternative 'quantization' the 
operator O has UV scaling dimension A_ . 
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Ion-shell- The Legendie transform must then be implemented to obtain the effective 
action F. In the alternative 'quantization', which is possible in the mass range 

e ~ t)' on-shell gravitational action is automatically a function of the 
vev (O) and expresses directly the effective action F. 

For static and homogeneous configurations the effective action F is proportional 
to the effective potential Vefj 

r[{0)] = -pV,.,V,ff{{0)) (3.5) 

where /3 is the period of the Wick rotated time direction and Vd-i is the volume 
of the d — 1 spatial directions. The vacua of the theory are easily determined by 
extremizing the effective potential, namely by solving the equations 

. (3.6) 



d{0) 



A vacuum is stable if and only if it is a minimum of the effective potential. 

In the rest of this section we use the language of the previous section to find 
general expressions of the on-shell gravitational action at finite temperature and 
charge density. 

3.2 Free energy from gravity 

The on-shell gravitational action / has a bulk contribution Ie and a boundary con- 
tribution Igh from the Gibbons-Hawking term. Here we evaluate with a UV cutoff 
ai u = Uq. a standard computation gives the free energy 

^ = Ion-shell = Mp^pVa-i (-W + /) e*^^" . (3.7) 

Implementing the equations (2.14) and (2.16) we obtain the following expressions, 
which will play a central role in this paper. 



or equivalently. 



^ = -e-oH^(<^„) _ 4vrT5 + / ' ^^J^, • (3.9) 



By taking the RG scale M = e"^° — oo we recover the free energy as a function of 
the source (in the case of the standard quantization) or directly the effective potential 
as a function of the vev (in the case of the alternative quantization). This step and 
the resulting expressions will be discussed in sections 5 and 6. 
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4. Renormalization group invariance 



Holography provides a direct way of accessing the concept of RG running, [6], and 
the analogue of /3-functions, [3]. In our case the /3 function for the scalar coupling 
can be obtained from (2.10), (2.11) to be 



d(/> 



dlogM dA 



m , m = -'^{d-i)d^\ogw{ct)) 



(4.1) 



where we identified A, in the standard manner, as the logarithm of the RG scale. 

The effective potential for the sources, evaluated at m = mq, is a function of the 
RG scale M = and the sources 0o, but does not depend explicitly on uq, namely 



duo 



. 



(4.2) 



The RG running of the effective potential is controlled by the /3-functions (4.1" 

d^ 



dAn 



dAo 



(4.3) 



which can be expressed further in terms of the superpotential function W. 
For instance, in the zero-temperature case 



dJ^ 



dAo d^c 



^dAo 



dW + W 



dAo 



dW - 2(rf - 1 



12 



w 



The superpotential equation (2.20) becomes at zero temperature and density 



dW -2{d-l 



W 



12 



w 



4(rf-i: 



V 

w 



Consequently, 



-^ = -4(ci-l)e'^^o^ 
dAo ^ ' W 



(4.4) 



(4.5) 



(4.6) 



5. Vanishing charge density 

The case of zero temperature and density is well studied. Expressions for the effective 
action in this case can be found for example in [7, 10, 11]. By setting T = p = in 
the general equation (3.9) we find 



(5.1) 



which expresses the free energy directly in terms of the 'superpotential' function W , 
[6, 12, 7]. In this section, we proceed to generalize this result at finite temperature 
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keeping the charge density p zero. At a finite cutoff Uq the general relation (3.9) 
becomes 

^ = _e'^^o|^(0^) _ 47rr5(0o) • (5.2) 

The temperature enters in this expression in two ways. First, it enters as a trivial 
overall factor in front of the second term in the rhs. Second, it enters non-trivially 
through (ph (the value of the scalar field on the horizon), which appears explicitly in 
the definition of S in eq. (2.17). 

Furthermore, one can derive the following useful identities between T, a, and 
the constant D. By differentiating and combining the equations of motion (specif- 
ically, the equations (5.13a), (5.13b) and (5.13c)) one can express the function / 
as 



2V + WDe-'^^ 



2{d-~l) 

Consequently, at the horizon we find that the following equation must hold 

2Vh + WhDe-''^^ = , (5.4) 

where we have defined 

Vh := V^(0/.) , := W{(Ph) , A := A(0O • (5-5) 
Now we can use eq. (5.4) together with the relation 

D = -AttTc'^^-^^^'^ (5.6) 

to eliminate e^'^ 

gA, ^ 27rT^ (5.7) 
Vh 

and express D in the computationally more convenient form 

D = -2{2nTY . (5.8) 
This allows us to recast eq. (3.8) (or (5.2)) into the form 

7 = -e'^°Wi<f)o) - 2{2'KTf f ^) . (5.9) 



0- 



The above equations express the free energy at a finite cutoff uq where : 
In order to remove the cutoff and make contact with the boundary QFT data, J and 
{O), we should scale Uq — )■ —oo and Aq — )■ +00. At the asymptotic boundary, where 
M — > — 00, the leading order terms in the expansion of the scalar field are 

= ae"'^- + ... + Ce"^+ + ... . (5.10) 
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Then, in the case of the standard quantization, the holographic dictionary dictates 
the relation 

a = J, C={0) (5.11) 

and, according to the general relations (3.2), the asymptotic coefficients a, ( obey 
the equations 

( = - — or equivalently a = . (5-12) 

da dC 

This relation together with a UV boundary condition, that provides an additional 
equation between a and (, determines the vacuum of the theory (and the full bulk 
solution) completely. 

In the case of the alternative quantization the roles of leading and subleading 
coefficients is exchanged. We proceed assuming the standard quantization. 

5.1 UV expansions 

In the case at hand, a solution of the bulk equations is determined by solving the 
system of ffist order equations^ 



W") f- WW'f = 2V , (5.13a) 



2{d-l] 

W'f = De-'^^ , (5.13b) 
1 ~W 

Useful information can be obtained by solving these equations perturbatively 
near = in the UV AdS asymptotic region. Given the expansion (2.2) we set 



oo 



n=0 n=l n=l 

where Wn, fn,gn have a separate expansion of the form 



oo 



Wn, fn, 9n = J2 • (^•^^) 



m=0 



Later we will fix 5 = We have assumed A_ < |. In the special case where 
A_ = A+ = I the latter expansion is performed in powers of log instead of 0. We 
will not consider explicitly this case here (see, however, [14] for a related discussion). 



^There is another class of potentials with exponential asymptotics related to decompactification 
limits in string theory, [1] . The superpotential for such asymptotics was analyzed in detail at finite 
temperature in [f3] and at finite temperature and density in [1]. 
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Up to next-to-leading order we obtain 



dW, 



TJ/'2 



- Wo = 2V , (5.16a) 



2{d-l) 

W;,f[ = De-'^'^g, , (5.16b) 
j^WqWi - 2W^W[ + 2Vfi - WqW'J[ = . (5.16c) 

The solution to the first equation, Wo, is the same as W in the zero temperature 
case. Soon we will see that Wq is the term that controls the UV divergences. In that 
respect, (5.16a) is in agreement with the fact that UV divergences are insensitive to 
the temperature. It is known that there are two separate solutions to (5.16a) with 
perturbative expansion 

wt' = - 1) + + f :r4^ - T,^7:r^^-7^) + . (5.17) 



2 ^ \d-4:A± 16(rf- l)(rf-4A±) 

It is also known [7, 10] that the Wq^"^ solutions do not allow for non-zero sources, 
whereas the H^o"'' ones do. Hence, in what follows we will restrict our discussion to 
the Wq~^ solutions and will drop the superscript 
The next order correction Wi takes the form 



m^wr"'-H«.y^^d,^-^-Be-^J. (5.18, 

The first term on the rhs of this equation, W^ , is the zero-temperature value of 
Wi 

„^ d r / rlfrl - 9A ^ Orliu \ _ .1 

(5.19) 



^^^lii^^^m^r^^^'' ''''''' 



and the second term a temperature-dependent contribution whose explicit form is 
determined in appendix B.l. In (5.19) C is a constant fixed by infrared regularity to 
a particular model-dependent value that we will denote as C*. 

5.2 The UV region Aq — )• +00 and renormalization 

As uq — )■ —00, 00 behaves at leading order as (see eq. (5.10)) 

00 = ae"«'^- + ... . (5.20) 

At the same time, in order to have a regular limit for 

1 1 



= e^M ^ '{! + ...) = e^»+"« [^] {! + ...) (5.21 
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we require Aq = —uq — )■ +00. 

Then, one can check the following limits 

lim e'^^Wo = 00 , (5.22a) 
lim e-^^Wi = C,a^ - '^—r^D , (5.22b) 

uo-)--oo d 

lim e'^^W„>i = . (5.22c) 

The divergence in the free energy arises solely from the first term (5.22a) that can be 
removed by subtracting the zero-temperature superpotential, W^c|j._q! for arbitrary 
constant C. This subtraction removes the divergent piece and shifts the constant C 

d 

of the a^- term in Wi. Therefore the renormalized superpotential is 

iy«(0, T) = W{C,, T, 0) - Wo{C, 0, 0) (5.23) 

Note that the renormalized superpotential Wr depends on the arbitrary parameter C 
used in the subtraction. This is a standard "scheme dependence" as in renormalized 
QFT. 

Using the expansion (B.2) we find that as we remove the cutoff the function 

_i r't'h w 

S = e 2 •''^o takes the form 

S = #0:^e^(^^) (5.24) 
where V{(j)) is a function with an analytic expansion in powers of 0. Then, 

and eq. (5.4) becomes 



27rTa ^ = 0, e^^^^^")-^ . (5.26) 

Wh 

The significance of this relation is the following. By inverting it we establish that (ph 

]_ 

is a function of the dimensionally proper combination 27rTa . 

We conclude that when we send Aq — > +00 (and subtract the divergence) the 
free energy (5.9) at finite temperature (expressed in terms of bare UV variables) 
takes the form 

^(a)^(C-C.)a--^(2.r)^(f)"'. (5.27) 

Hence, in the standard quantization the effective potential K//(C) can be determined 
from the Legendre transform 

VeffiO = -{C - a)a^ + {2nTY (^Y \ aC , (5.28a) 
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C=^. (6.28b) 

In the alternative quantization a represents the vev of the dual operator and the rhs 
of eq. (5.27) expresses directly the effective potential. These expressions are one of 
the main results of this paper. 

d 

Notice that the T = part of the free energy (5.27), proportional to a^- , follows 
also from the scale invariance of the planar boost-invariant solution and agrees with 
the expressions in [11]. We will discuss the temperature dependence further in section 
7. 



6. Finite temperature at finite density 

In this section we extend the discussion to arbitrary temperature and charge density. 

We are now solving the more complicated system of first and second order differ- 
ential equations (2.21a)-(2.21c). The expansions (5.14) have a slightly more involved 
form which is presented in appendix B.2 along with the corresponding expansion of 
the equations of motion. 

One can show that the analog of eq. (5.3) is 

^ _ 2^ - + ^-^^^ + It 

2(d-l) 

Hence, at the horizon j{uh) = implies 

2^'. - ^5(5^ + ^-"'W'. + £ ^Ur^) = . (6.2) 
Implementing eq. (2.16) we obtain 

2K - - ^^TSe-'^^W, = . (6.3) 

Finally, with the use of eqs. (2.17), (5.20), (5.24) we find 

1 2(d-l) ^(''-^) 1 -J- p(0 ) 

2Vh - —p^a ^- 4>^ - e-2^('^'^) - iTiTa ^- e-^Wh = . (6.4) 

V{(f)) is not identical to the function defined in eq. (5.24) for p = 0, but is defined in 
the same way. Equation (6.4) is the generalization of eq. (5.26). We conclude that 
4>h is now a function of the combinations 

n_ ^ 2{d-i) 

Ta ^- , p^a . (6.5) 
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The free energy is given by eq. (3.9) 

As we remove the cutoff by taking Aq — +oo, and appropriately regulate the diver- 
gences by subtraction as in the zero-charge case, we obtain 

- lim W{(j)o)e'''''> = {C- a)a^ + . (6.7) 

Ao-5>oo a 

The sole contribution comes from Wi (see eq. (B.9a)) as in the p = case. Recall, 
however, that D is now given in terms of the temperature by eq. (2.16). Hence, we 
obtain 



^ p 9W _ 1 9/7 _ 1 /■'^o 1 



2rf - 1 . „ _ 2d - 1 o Z"^" 
Using eq. (5.24) and defining for convenience the function 



'"00 ]^ 

•^•('^^) - / :rawTi77 (6-9) 



^id-2)AZW' 

h 



we find the expression 



a a 

The functions P, are in general complicated model-dependent functions that can 
be determined with the use of numerical methods. 



7. Scaling asymptotics and phase transitions 

The unwieldy model-dependent behavior of the general expression (6.10) simplifies 
and reveals universal information in the vicinity of critical points. In this section 
we discuss transitions that involve the 'normal' phase with vanishing scalar conden- 
sate. For simplicity, we consider the case of the alternative quantization in which a 
captures the vev of the dual operator and the equations (5.27), (6.10) express V^ff 
directly clS cl function of a. 

7.1 Vanishing charge density 

As an interesting warmup exercise we first consider the case of vanishing charge 
density. 

Solutions with everywhere small values of the scalar field (f) are big black holes 
with horizon in the UV region, namely 0/j ^ 1. For such solutions the UV perturba- 
tive expansions of subsection 5.1 are valid for the whole solution outside the horizon 
and one can use them to evaluate perturbatively the rhs of eq. (5.27). 
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First, inverting the relation (5.26) we find at leading order 

(if)'^..,... 



1 



The dots indicate subleading terms in inverse powers of Ta ^- . From this expression 
we learn that in general 0/i ^ 1 can be interpreted either as the small vev limit at 
finite temperature or as the high temperature limit at fixed vev. 

Inserting (7.1) into the expression for the effective potential and expanding up 
to quadratic order we find 

V;„ = -Ca--(2d-l)(— ) + 

(7^2) 

where, compared to previous formulae, we have replaced the constant C — )■ C+C^:. In 
the second, (9(q!°), term we recognize the standard T'^ free energy of a d-dimensional 
conformal field theory. 

In the presence of a double-trace deformation on the field theory side 

6C ~ gO^ (7.3) 

the effective potential at zero temperature becomes [11] 



Veffia)L.=ga'-Ca^- . (7.4) 



Hence, assuming for concreteness C < 0, we learn that the normal vacuum at a = 
is unstable when g < and a stable symmetry-breaking vacuum exists with vev 

2gA- 



"=(^j ■ P^^' 

Adding temperature in the presence of the double-trace deformation we obtain 
the effective potential 

K//(a) = -Ca^ - ET'' + g.fja^ + . . . (7.6) 

where gejf is the temperature-shifted effective double-trace coupling 

g,ff = g + GT''-^''- (7.7) 

and -E, G are positive (a, T)-independent constants that can be read off eq. (7.2). 

The stability of the normal state at a = can be determined immediately from 
the sign of geff- In particular, the normal vacuum becomes unstable when gefj < 0. 
The critical temperature that separates the stable from the unstable regime is 
obtained by setting 

.e// = ^ T,= f-|,)^ . (7.8) 



a, 

This formula is in agreement with the result obtained in a different manner in [8]. 
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7.2 Finite temperature at finite density 

A similar analysis of the stability of the normal phase can be performed in the 
more general case of finite temperature and density using the formulae of section 
6. Expanding the effective potential up to quadratic order 0{a'^) (in the presence 
of a double-trace deformation) we obtain a formula analogous to (7.6) with a more 
complicated effective double-trace coupling. The explicit formula appears in eq. (C.9) 
of appendix C. 

A finite-temperature transition at fixed p can be determined again by setting 
Qeff = 0. This transition occurs at a quantum critical point when the equation 
Qeff = can be solved for T = 0. This is possible at the critical double-trace 
coupling 



whose explicit p-dependence, and the constants Ci,C2 are determined in appendix 
C. To the best of our knowledge this analytic expression is new. In the vicinity of 
the quantum critical point we observe the following scaling of the vev 



8. Outlook 

In quantum field theory the effective action contains all the information needed to 
determine the vacuum structure and generic correlation functions of the theory. This 
information can be used, for example, to identify critical points and phase transi- 
tions at finite temperature and/or finite density, and as such it is very useful for 
descriptions of real world systems in high energy physics or condensed matter. Un- 
fortunately, the direct computation of the full effective action of interacting QFTs 
is typically an almost impossible task. As a result, the available effective theories 
are limited to descriptions of the vicinity of special points in parameter space. For 
instance, once a critical point and an order parameter that characterizes the corre- 
sponding transition is identified, a Landau- Ginzburg description can be employed to 
describe the properties of this order parameter in the vicinity of the critical point. 

Theories with a weakly curved dual gravitational description allow us to go be- 
yond these limitations and derive (with specific rules) the effective action of an order 
parameter in generic regimes of parameters. In this work we have initiated a sys- 
tematic study of the effective action of theories with a dual description that involves 
a class of Einstein-Maxwell-dilaton theories. We focused on spacetime-homogeneous 
configurations and the corresponding effective potential at finite temperature and 
charge density. The resulting formulae extend the usual Landau-Ginzburg treat- 
ments away from critical points and allow us to probe the full off-critical behavior 




(7.9) 



A 



{O) ~ {g, - g) 



(7.10) 
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of the system. The latter is a rather comphcated model-dependent function of the 
parameters. As expected, the behavior simplifies near critical points where we can 
recover specific critical exponents and other data of the transition. 

There are several interesting extensions of this work that are worth exploring 
further. First, it would be useful to explore the relation of our approach with the 
more systematic Hamilton- Jacobi study of holographic renormalization [6, 7]. For 
example, a more thorough understanding of the precise relation between the super- 
potential function W that we introduced and Hamilton's principal function at finite 
temperature and density would be useful. 

Second, for many applications it is useful to know the full effective action in- 
cluding derivative interactions, a subject interesting for applications both to finite 
density systems and cosmology. Work in this direction using the Hamilton- Jacobi 
formalism has appeared in [10, 7], where explicit expressions at zero-temperature and 
density are provided. In this work we focused on effective potentials for scalar order 
parameters, but more generally effective actions for generic tensor order parameters 
{e.g. vectors) are of intrinsic interest in various applications. 

The Einstein-Maxwell-dilaton theories in this paper have appeared in numerous 
applications of the holographic duality to QCD and condensed matter theory. A 
class of these applications in condensed matter theory refers to holographic super- 
conductors [8] and holographic models of magnetism [15]. The effective potential 
that we analyzed above can be used to study the vacuum structure of these, and 
more general, systems. These systems can also be used as basic building blocks in 
the construction of corresponding Josephson junction networks within the framework 
of designer multi-gravity theories [16]. Applying the results obtained above to this 
context we can determine the effective potential of corresponding network theories 
and use it to study their finite temperature and density physics. 
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Appendices 

A. The Einstein-abelian-Higgs model 

As explained in section 2.2 the action of the Einstein-abelian-Higgs model (2.24) 
reduces to the action (2.25) 



d-l 



xy/-g 



R 



Z{<t>) ,2 



-(d<pf-q^J{<P)A' + V{<P) 



+ Igh (A.l) 
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by setting 9^6' = 0. (/> is a real scalar field that represents the modulus of the complex 
scalar field with U{1) charge q of the Einstein-abelian-Higgs model. In the present 
appendix we summarize the equations of motion of the action A.l and the form they 
take when a superpotential function W is introduced as in section 2.1. 
In the domain wall frame 

d?/2 

ds^ = e2^(") (-f(u)dt^ + dxMx,) + -— , A = At(u)dt , = (A.2) 

f{u) 

we obtain the following independent set of second order differential equations 

d / ^id-2)A2^\ _ 2q^e^'i-^)^-Ut = , (A.3a) 
du \ / / 

2(.-lM + ,^ + 2,V-^.-0. (A.3b) 

{d - l)Af + (^d{d - l)A^ - ^02^ f-V+ ^Ze"2^i2 _ ^2 jg-2A^-1^2 ^ Q ^ (A 

We introduce the function W by requiring the flow equation 

<j) = W' . (A.4) 

Then, eq. (A. 3b) gives 

where we have chosen a particular integration constant that fixes the definition of 
W. Eq. (A. 3c) becomes 

W'iW'fy - (W + 2q^ J' "j^A^^ W'f - e-^^ZA\ - 2q\--^^-A\ = 

(A.6) 



and eq. (A. 3d) 



2 

-Ze-^^i? _ 2gVe-2^/-^A2 = . (A.7) 
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Eq. (A. 6) can be integrated to obtain the expression 



W'f 



D + e'^^H 



-2A 



d0 



HW 



(A. 



where H{(j)) is the function that solves the differential equation 



H' 



H 



2{d-l) W 
and D is an integration constant. Using the horizon regularity condition 



AttT 



(A.9) 



(A.IO) 



U=Uh 



we determine D as follows 



</>0 



-47rTe(P-^)^("'') = -^ixTS . (A.ll) 



The free energy takes the 'universal' form (3.8) 



T 



(A.12) 



B. UV expansions 



B.l UV expansions for vanishing charge density 

In this appendix we consider in more detail the solution of the equations (5.16b) 
(5.16c) 



-^W^W^ - 2W;,W[ + 2Vh - W,W'J[ = 
From eq. (5.13c) and the expansion 
d 



(B.la) 
(B.lb) 



2{d 



d 



1 + 



2(rf-l)4 n 

1 , 92 Qs gl . 
log^i + — + + • • • 

9\ 9i m 



/2 



+ 



we obtain 



5-1 = e 



<k r'i' 

2{d-l) J0o ^w/i 



(B.2) 



(B.3a) 
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d 



.,W, (Wi_W[ 



etc. 



(B.3b) 



Wq is a solution of eq. (5.16a) with the (— ) expansion in (5.17). 
We can integrate eq. (5.16b) to find /i 



J4>o 



(B.4) 



Then eq. (5.16c) becomes 



d 



W[ -Wi = — fi -f{ 



with solution 



(B.5) 



De 



-dAo 



Wq 

2W^ 



(B.6) 



where W^^ is the zero-temperature value of Wi 



(T=0) 



1 + 



d{d - 2A. 



2dvA 



4(d- l)(ci-4A„) Alid-AA^ 



+ W 



(B.7) 



As was commented in the main text, C is a constant fixed by infrared regularity to 
a particular model-dependent value that we denote as C^,. 

B.2 UV expansions at finite charge density 

At finite charge density we are solving the more complicated system of first and 
second order differential equations (2.21a)-(2.21c). The expansions (5.14) are slightly 
more involved 



■ dn+2m(d-l) ^ 

n.m 

For the first few orders 
W = W0+W1+W1 + . . . , Wi = (f) 



ioT h= {W,f,n) . (B.8) 



2(d-l) 



Wi,nr , H^i = ^- J] W^i,n0" , (B.9a) 



n=0 



n=0 



7^ = 7^o+7^l+7^l+... , 7^l 
/ = /o + /i + 7i + ... , /i 



^7^l,„0^ 7^l = ^- ^7^l,„0^ (B.9b) 

n=0 n=0 

/l = 5^/l,n0V (B.9C) 



n=0 



n=0 
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The form of the equations and the corresponding solution at p = suggests setting 

7^o = 7^l = , /o = l . (B.IO) 

Inserting these expansions into the equations of motion (2.21a)-(2.21c) we find 
that the functions Wq, Wi, fi are the same as in the p = case and 

^, = e"^^^"^^"«<^, (B.lla) 



where gi is the function that appears in eq. (5.14) and 0* is an integration constant. 



C. Small-0/i expansion details 

In the regime 0/i ^ 1 we have the following expansion 

(l)h = Aia{l + A2a^ + . . .) , (C.l) 
V{<j)H) = B,ct>l + . . . , (C.2) 
F.{<Ph) = (1 + C24>1 + ■■) (C.3) 



Bi can be determined from the definition (5.24) 

- "32(rf-4A_) • ^^-^^ 

The coefficients Ai, are determined by expanding eq. (6.4) in (^h- More specif- 
ically, Ai is determined by solving the algebraic equation 

2(d-l) 1 

p%^- + 87r(d - l)TAf- - 2d(d - 1) = . (C.5) 
is determined in terms of Ai as follows 

(2rf(rf - 1)^1 + m2)A2 + 25ip2A/^^^ - 27rTAf^^^(4(c/ - l)^i + (A_ - AB^)) 

^ 2(d-l) ^~ • 

p2Ai^- +47rTAf- 
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Finally, the coefficients Ci, C2 are determined by expanding the expression (6.9) 
and taking the 0o — > limit. We find 



Co 



1 d - 16A. 



1 



d-2 
3 



Vi - 



(C.7) 



(C. 



16rf-4A_ ' A_(rf-4A_) V 16(rf-l), 

Expanding the effective potential up to quadratic order O(a^) (in the presence 
of a double-trace deformation) we obtain a formula analogous to (7.6) with a more 
complicated effective double-trace coupling 



9eff = 9 

( 



2d-\ 
d 



2{d-l) 



M{d-1) - p^A-^ 



2(d-l) 



=^ -{d- 1) ^ ^ 



v 



2d{d - 1) - 



2(d-l) 
A_ 



/ 



— / rl — 9 



A, 



(C.9) 



D. Non-linear Einstein-scalar actions 



As a qualitatively different example, in this appendix we consider actions of the 
general form 

I = Mp^ j d^^^x ^g[R + V{<P)F {{d<pf)] + Igh ■ (D.l) 

is a real scalar field. For the function F{x) we demand -F(O) = 1, -F'(O) = — y < 
so that in the long-wavelength limit we obtain the canonical scalar field action 



R + v{x)--{dxf 



X = U dx^/V{x) , V{x) := Vi<Pix)) ■ 
Jo 

A commonly encountered example is the DBI case with 

Fdbi{x) = y/l-x . 



(D-2) 



(D.3) 



In what follows we describe how the equations of motion of this system can be 
reduced to a first order system for zero-temperature solutions. A general treatment 
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of the finite temperature case lies outside tlie immediate scope of this paper. It 
can be performed along the lines of the previous subsection. Another interesting 
generalization suggested by the DBI case would be to include a U{1) gauge field. 
Working again in the domain wall frame we set 

ds^ = e2^(")(-dt2 + dx'dx') + du\ (j) = ^{u) . (D.4) 

The equations of motion reduce to the following set of independent differential equa- 
tions for the unknown functions A{u), (f){u) 

{d - 1)A = , (D.5a) 

d{d - l)A^ -VF + = o . (D.5b) 

y is a function of 0, F a function of 0^ and we are using the notation 

In this case the superpotential ansatz (2.10), (2.11) is generalized to 

Inserting this ansatz into the equations of motion (D.5a), (D.5b) we obtain the first 
order system 

^ _ _ _ , (D.8a) 



4(d- 1 

W = -2VFW . (D.8b) 

As an illustrating example consider the DBI case (D.3). One can solve explicitly 
the second equation (D.8b) to obtain 

~ eW V 

W = , , FnBi = , , e = ±1 . (D.9) 

Substituting into the first equation (D.Sa) we find a single non-linear first order 
equation for W 

eW^ + V^= ,f , W^Jv^ + W^. (D.IO) 
4(a — 1) ^ 

One of the novelties of the non-linear case is the presence of multiple branches of 
solutions parametrized by the free parameter e = ±1. 

The same results are obtained easily with the use of the Hamilton- Jacobi for- 
malism. For static and homogeneous configurations Hamilton's function is 

Sh = -Mp^ I d'^x^/^W{(t)) = -M^-^/3\/d„ie'^^W(0o) . (D.ll) 

J u=un 
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Eq. (D.8a) is the Hamilton- Jacobi equation and eq. (D.8b) is the first order flow 
equation for the scalar fleld 0. 

The formula (D.ll) in the Hamilton- Jacobi formalism implies that the (rescaled) 
free energy is 



a result that can also be verifled easily by direct computation. 
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